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ABSTRACT 

The  singularity  expansion  method  (SEM)  has  been  applied  to 
determine  natural  resonances  of  a  perpendicular  crossed  wire  structure 
parallel  to  an  Imperfect  ground  plane.  To  account  for  the  finite  con¬ 
ductivity  of  the  ground,  the  Fresnel  reflection  coefficient  method  has 
been  used  to  follow  the  trajectories  of  the  natural  resonances  as  the 
conductivity  and  relative  permittivity  of  the  ground  vary. 


I.  INTRODUCTION 


The  general  problem  of  the  electromagnetic  Interaction  of  wire 
structures  and  external  sources  of  excitation  has  been  of  Interest  for 
many  years  [1,2].  In  a  series  of  studies  completed  In  recent  years  the 
progression  has  beqa  as  follows:  the  single  isolated  cylinder  in  free 
space  in  the  frequ*;icy  domain  [3].  a  perpendicular  crossed  wire  struc¬ 
ture  in  free  space  in  the  frequency  domain  [4],  a  non-perpendicular 
set  of  crossed  wires  in  free  space  in  the  frequency  domain  [5],  a  study 
of  a  cylinder  [6]  and  a  set  of  perpendicular  crossed  wires  [7]  using  the 
singularity  expansion  method  (SEM)  [8],  a  cylinder  [9]  and  a  set  of  per¬ 
pendicular  crossed  wires  [10]  over  a  perfectly  conducting  ground  plane 
using  image  theory  and  SEM.  Recently  two  studies  of  a  cylindrical 
scatterer  have  considered  different  effects  of  the  ground:  in  the  first 
the  effects  of  a  close  proximity  perfectly  conducting  ground  plane  on 
the  clrcimferential  distribution  of  currents  [11]  have  been  studied,  and 
in  the  second  the  effects  of  a  finite  conducting  ground  have  been  studied 
using  SEM  and  Fresnel  reflection  coefficients  [12].  Finally  preliminary 
studies  were  performed  to  analyze  the  crossed  wire  and  the  cylinder  over 
the  finite  conducting  ground  using  the  Sommerfeld  formulation  [13,  14,  15] 
in  the  frequency  domain  [16].  Simultaneously,  others  were  applying  trans¬ 
mission  theory  to  gain  information  on  this  same  problem  [17]. 

The  problem  considered  in  this  study  is  that  of  determining  the 
behavior  of  the  natural  resonances  in  SEM  for  a  set  of  perpendicular 
crossed  wires  over  an  Imperfect  ground  using  the  Fresnel  reflection 
coefficient  method.  While  the  technique  has  inherent  limitations,  trends 


can  be  established  that  are  of  value  in  themselves  as  well  as  providing 
guidelines  for  more  correct  formulations;  e.g. >  e  Soamerfeld  treatment. 
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II.  FORMULATION 

This  analysis  deals  with  Che  Interaction  of  a  set  of  perpendicular 
crossed  wires  whose  plane  is  parallel  to  an  imperfect  ground.  The  effect 
of  the  ground  on  the  wire  structure  is  approximated  by  the  Fresnel  reflec¬ 
tion  coefficient  method.  SEM  is  used  to  follow  the  trajectories  of  the 
natural  resonances  as  the  ground  conductivity  (Og)  and  the  relative 
ground  permittivity  (cr)  vary.  In  analyzing  the  scattering  problem,  thin 
wire  approximations  are  utilized.  That  is,  the  wires  are  required  to 
satisfy  the  follow: ng: 

(1)  the  current  and  charge  densities  are  constant  around  the 
periphery  of  the  wire  cross  section, 

(2)  the  axial  currents  are  zero  at  the  free  ends  of  the  wires, 

(3)  ka  «  1, 

(4)  a  «  i 

where  k  is  the  propagation  constant  defined  by  2w /wave length,  a  is  the 
radius  of  the  wire,  and  1  is  the  length  of  the  wire. 

Graves  [17]  has  applied  Pocklington  theory  to  a  set  of  orthogonal 
crossed  wires.  Figure  2-1.  If,  in  the  method  of  moments  solution  to  the 
thin  wire  Pocklington  type  integral  equations,  sinusoidal  pulse  expansions 
are  used  for  the  unknown  current  distributions  and  delta  testing  functions 

J 

are  used,  relatively  simple  results  are  obtained  that  are  equivalent  to 
the  following.  If  the  current  distribution  on  a  filament  is  assumed  sinus¬ 
oidal,  exact  analytical  expressions  for  the  accompanying  field  components 
are  [18] 

r2l  •  ^(t)  -  |  !’(?!>  Cff.fj)  -  I'(*2>  G(t,*2)j  (2.1) 

for  the  axial  electric  field  where  s  is  the  complex  Laplacian  frequency, 
s  ■  a  ♦  jo);  the  current  filament  extends  from  to  ^2,  I'  the  axial 


v.  eg 


derivative  of  the  current  and 
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*2  - 

Ift-ni 


G(r,?i) 


e-s|r-ri|/c 


The  corresponding  radial  component  of  the  electric  field  is 


;2i  •  *cn 

-«  <?2  > 


1 

4ttesp2i 


(r-?2)I'(r2)  G(r,r2) 


( |r-r2 | )G(r,r2)  -^(rj) ( | r-rj | )G(r 


*A>]  | 


(2.2) 


and  p2^  is  the  radial  variable  or  the  perpendicular  distance  from  the  axis 
of  the  current  filament  to  the  field  point. 

To  derive  the  system  of  equations  for  the  current  distribution  on  the 
crossed  wires  by  the  method  of  moments  the  various  wires  are  divided  into 
segments.  The  current  on  the  kth  segment  of  the  mtl*  wire  is 


Wx> 


a^sinh 


[s(Xlc+i  -  x)/c] 


+  °jt+^sinh 


^s(x-xk)/cj 


(2.3) 


8inh[,(Vl  -  ^)/c] 


It  follows  that  is  the  axial  current  at  the  point  x  ■  x^  where  x  is  a 
general  variable  along  any  given  wire. 

For  crossed  wires  in  free  space  the  procedure  is  as  follows.  Assume 
currents  Iiv(y)  and  I2k(x)  flow  on  the  wires  in  Figure  2-1.  Substituting 
(2.3)  into  (2.1)  and  (2.2)  leads  to  electric  field  expressions  that  can  be 
evaluated  at  points  on  the  conducting  structure.  For  example,  a  current 
zone  on  wire  1  will  induce  an  axial  field  at  other  zones  on  wire  1  through 
(2.1)  and  an  axial  field  (as  far  as  wire  2  is  concerned)  at  zones  on  wire 
2  through  (2.2).  The  total  axial  field  at  each  element  due  to  the  induced 
currents  is  set  equal  to  (-^incident ) tangent  at  a  discrete  set  of  match 


points  along  each  wire,  and  solutions  for  the  a's  may  be  obtained.  In  the 
remainder  of  this  report  the  field  contributions  (2.1)  and  (2.2)  are  called 
primary  scattered  (p.s.)  excitation  terms. 

The  effects  of  a  perfect  ground  plane  may  be  accounted  for  by  introducing 
image  currents  (parallel  to  the  ground  plane  and  opposite  in  direction  to 
the  object  currents)  and  calculating  their  E-field  contributions  at  match 
points  on  the  crossed  wire  structure.  A  modification  on  the  G  functions  in 
(2.1)  and  (2.2)  is  necessary  due  to  the  distance  factor  now  depending  on  the 
height  above  the  ground,  but  the  procedure  is  straight  forward  [10].  In  the 
perfect  ground  plane  problem  the  tangential  component  of  the  incident  field 
is  now  added  to  the  two  p.s.  terms  plus  two  image  terms  to  produce  a  total 
tangential  E  field  on  the  wires  that  equals  zero.  The  image  terms  are 
referred  to  as  secondary  scattered  (s.s.)  terms. 

Riggs  and  Shumpert  [12]  have  recently  studied  the  problem  of  a  cylinder 
over  an  imperfect  ground  by  modifying  the  image  term.  The  modification  is 
a  multiplicative  factor  on  the  image  term  that  equals  the  Fresnel  reflection 
coefficient  for  a  non-magnetic,  finite  conducting  ground  with  conductivity 
0^  and  relative  permittivity  Limitations  do  exist  in  this  approach: 

the  Fresnel  reflection  coefficients  are  derived  for  incident  plane  waves, 
thus  the  "scattered"  field  at  the  interface  between  the  lossy  ground  and 
free  space  should  approximate  a  plane  wave.  Sarkar  and  Strait  [15]  have 
studied  this  method  in  the  frequency  domain  and  observed  that  their  results 
are  within  10%  of  the  Sommerfeld  formulation  if  the  dipole  was  at  least 
(0.25  W/£r)  above  the  ground. 

Consider  the  geometry  of  Figure  2-2.  The  current  at  82'  contributes 
to  the  scattered  field  at  S2  due  to  the  p.s.  term  and  due  to  the  s.s.  term. 
Since  the  discussion  is  limited  to  axial  currents  only  the  axial  field 
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component  from  82’  will  contribute  to  the  p.s.  and  s.s.  terms.  It  follows 
therefore  that  s.s.  field  lies  in  the  plane  of  incidence,  and  the  Fresnel 
reflection  coefficient  becomes  (19) 


RPA  - 


(£r  +  X)  sirajj (a; ,sj)  -  V (cr  +  X)-cos2ij/(s?  .sj) 

(er  +  X)  siroKs2,82)  +  "V^er  +  X)  _  cos2  i(>(s2,  sp' 


(2.4) 


where  X  ■  120  iTOg  ,  y  - 

Y 

The  particular  sign  associated  with  RPA  depends  on  the  choice  of  the  direction 
of  the  reflected  field  (Riggs  and  Shumpert  use  one  convention  and  this  study 
uses  a  different  one,  hence  our  reflection  coefficient  differs  by  a  multi¬ 
plicative  factor,  -1).  In  the  limit  as  Og-*-  ”,  image  theory  and  the  reflection 
method  become  the  same. 

The  crossed  wire  problem  is  geometrically  more  complex  but  can  be 
treated  by  the  reflection  method.  Consider  the  current  element  at  sj;  see 
Figure  2-3.  The  ground  scattered  field  is  neither  parallel  nor  perpendicular 
to  the  incident  plane.  Thus  the  axial  field  Eg^(sp  is  resolved  into  two 
terms,  Esj(sp()  and  Esl(si^J_*  that  are  parallel  and  perpendicular  to  the 
incident  plane,  respectively.  Similarly,  the  radial  field  component, 

Ep(sJ),  is  resolved  intp  parallel  and  perpendicular  terms.  First,  Ep(s|) 
is  projected  into  the  Sj.s^  plane  or 

Ep2*  1  “  Ep(s[)  sin  6  » 

see  Figures  2-3  and  2-4.  Then  the  components  of  EpS2’sl(s[)  perpendicular 
and  parallel  to  the  incident  plane  are  found.  Figure  2-4.  Finally,  the 
components  of  Es j (sj)  are  defined;  see  Figure  2-5. 


a 
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The  total  scattered  field  in  the  plane  of  incidence  at  the  scattering 
point  is 

Ep(sP||.  "  Esi(8P|I 

and  the  total  reflected  "parallel"  field  at  s^  is 

RPAlEP<8i>||  -  Esi<8i>||]  (2-5) 

Similarly  the  total  reflected  "perpendicular"  field  at  s^  is 

Rve[ep(s'>]_  +  E8l(apl]  (2.6) 

In  (2.S)  and  (2.6) 


(2.6) 


RpA  _  (er  +  X)  siniKs2.si)  -  V(£r+  X)  -  cos2^(a2»8{) 

(er+  X)  sin4i(82»sp  +  X)  ••  co82<Ka2»8p* 


m  8in!(32,8[)  -V(cr  +  X)  -  co824»(82,sl) 


sini(32»s[)  +V(£r  +  X)  -  co82iKs2,sp 


(2.7) 


(2.8) 


The  total  scattered  field  at  a  point  on  wire  2»  e.g.  Sj,  can  now  be 


written  as  the  vector  sub  of 

1)  Eprimary<82’8P 

2)  Eprimary(82’8P 

3)  ^secondary (s2* 8 P 
^  Secondary  ^s2* 8  P 


(see  Figure ; 2—2) 


(see  Figure  2-3) 


(see  Figure  2-2) 


(see  Figure  2-3) 


Item  1  and  2  are  simply  the  free  space  terms.  Item  3  is  the  secondary 
scattered  term  of  Riggs  and  Shunpert  while  item  4  is  the  combination  of 
(2.5)  and  (2.6).  The  tangential  component  of  (this  resultant  vector  plus 


the  incident  field)  is  set  to  zero.  Thus,  a  system  of  linear  equations 


< 


13 


for  Che  unknown  current  expansion  coefficients  ie  produced  in  the  for* 

Z(s)  7(e)  -  E(a)  (2.9) 

where  Z(s)  is  the  square  matrix  of  the  system,  J(s)  is  the  coltsnn  matrix 
of  the  a's(see  2.3),  and  E(a)  la  a  colunn  matrix  whose  elements  are  the 
particular  values  of  the  Incident  field  components. 

Since  the  techniques  of  SEM  are  well-documented  [6,7,8,9,11,12],  the 
details  will  not  be  repeated.  The  problem  of  interest  la  the  behavior  of 
the  natural  resonances  as  the  ground  parameters,  Og  and  er,  vary.  To  find 
the  natural  resonances,  the  non-trivlal  solutions  to 

de-|z(a)|-0  (2.10) 

must  be  found.  Various  techniques  may  be  applied  (6»7,10);  In  this  study 
the  Muller  iteration  la  used. 
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III.  RESULTS  AMD  CONCLUSIONS 


The  geometry  of  the  problea  is  defined  in  Figure  2-1.  In  particular, 
the  peraaetere  choaen  are 

h/L  -  0.2,  L/a  -  20.0,  iJ/lj  -  0.5,  +  tj  -  2i2  -  L. 

The  trajectorlee  of  the  firet  eix  natural  resonancea  are  determined  ae  the 
ground  conductivity  varies  frotc  120  mhos/m  to  0.001  ahos/m,  and  the  relative 
permittivity  of  the  ground  varies  from  1  to  35.  Table  3-1  and  Figure  3-1 
give  the  comparison  between  the  results  of  this  study  for  a  high  conductivity 
(120  mhos/m)  and  cr  equal  one  (1)  to  a  previous  study  of  the  crossed  wires 
over  a  perfectly  conducting  ground  (10].  The  observed  differences,  mainly 
in  the  real  parts,  can  be  attributed  to  two  factors:  1)  the  distance  factors 
on  the  image  terms  [10]  in  the  secondary  scattered  terms  are  defined  in  a 
slightly  different  way;  2)  the  conductivity  is  not  infinite.  Figures  3-2 
through  3-7  present  families  of  curves  exhibiting  the  results  of  this  study. 
The  point  x  on  each  figure  is  the  image  theory  [10]  location  of  the  resonance. 
Also  on  each  figure  is  the  location  of  the  free-space  resonance, .  Con¬ 
sider  Figure  3-2:  each  curve  on  this  figure  represents  the  trajectory  of 
the  sy.1,1  resonance  for  a  particular  choice  of  er(l ,5,10,20,35) .  The 
points  on  a  given  curve  tepresent  the  value  of  the  sy,l,l  resonance  on  the 
particular  curve  for  values  of  Og  (A  •  120,  B  -  10,  C  -  1,  D  -  0.1, 

E  *  0.01,  F  ■  0.001  mhos/meter). 

Several  conclusions  can  now  be  drawn,  a.  For  a  highly  conducting 
ground  ('  120  ahos/a)  the  resonance  location  is  relatively  independent  of 
the  cr  value  and,  in  fact,  approaches  the  image  theory  result  in  each  case, 
b.  For  an  er  ■  1.0,  as  the  ground  conductivity  approaches  aero  the  resonance 
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TABLE  3-1 

Complex  Natural  Resonances  for  Crossed  Wire  Structure 
(ti  +  li  -  L  -  2lit  L/a  -  20,  fci/ti  -  0.5,  h/L  -  0.2) 


resonance 

perfect  ground 

imperfect  ground 
°-120g£!r  er  -  1 

sy.1,1 

-0.0513  +  J  2 . 361 

-0.0759  +  j 2. 3420 

as,l 

-0.0898  +  J2.592 

-0.1246  +  J2.5619 

®y,i,2 

-0.1021  +  J3.769 

-0.1296  t  33.7450 

sy,2,2 

-0.3909  +  J5.740 

-0.4659  +  J5.7155 

sy,3,l 

-0.6691  +  j7.581 

-0.7816  +  J7.5489 

as, 2 

-0.6478  +  j8.034 

-0.7579  +  j8.003 
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location  approaches  the  free  space  value,  c.  For  large  er  values,  the 
trajectory  of  a  resonance  covers  less  of  the  complex  plane  as  o  varies. 

O 

The  curves  of  Riggs  and  Shumpert  exhibit  the  same  behavior  but  for  a  very 
different  shape  parameter  (L/a  ■  200  in  their  study). 


-«>  .V*.  •  • 
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